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AN EFFECTIVE METHOD TO COMPUTE CLOSURE ORDERING FOR 
NILPOTENT ORBITS OF ^-REPRESENTATIONS 

W.A. DE GRAAF, E.B. VINBERG, AND O.S. YAKIMOVA 

Abstract. We develop an algorithm for computing the closure of a given nilpotent Go- 
orbit in 01, where Qi and Go are coming from a Z or a Z/7TiZ-grading = ® 0i of a simple 
complex Lie algebra g. 



1. Introduction 

O 

^ ' One of the main tasks of mathematics is to describe certain objects up to a certain equiv- 

^ ■ alence relation. Often this relation is given by an algebraic group action. Then equivalence 

classes are orbits and orbit closures correspond to degenerations of our objects. Thus, de- 

O I scribing orbits of algebraic actions, as well as deciding whether one orbit lies in the closure 

of another, is an important and interesting problem. However, this is possible only in a 
very few cases. One of these instances is provided by the ^-groups introduced by the second 
author in the seventies, see [34], [35]. 

Let G be a connected reductive complex algebraic group and g = LieG its Lie algebra. 

T^ • Let ^ be a diagonalisable automorphism of q that either defines a Z or a Z/mZ-grading 

^ ■ = 00i) where the grading components Qi are the eigenspace of 9. Note that Qq = q^ 

is the subset of ^-stable points. Let Gq G G he a. connected algebraic subgroup such that 
Lie Go = So- If ^ extends to an automorphism of G, then Gq = (G^)°. The group Gq is 
reductive and its natural action on Qi is called a 6 -representation; the group Go, together 

^ . with its action on Qi, is called a 9 -group. 

\Q ', An orbit Goa: C Qi is said to be semisimple if it is closed, and nilpotent if its closure 



< 



o 



^^ I Gox contains 0. This is the case if and only if x is semisimple (respectively nilpotent) as 

an element of q. The elements of Qi inherit the Jordan decomposition x = s -\- n from q. 
Besides, Go-orbits Go(s + n) with the semisimple part s being fixed up to conjugation are 
classified by the nilpotent orbits of the 6'-group coming from the pair (fl^, ^|gj, where 0s C 
is the centraliser of s (and it is a reductive Lie algebra) and 6\q^ is the restriction of 6' to 0^. 
This indicates that nilpotent orbits are especially interesting . The 6'-groups have several 
^ . remarkable properties, one of them is that there are only finitely many nilpotent Go-orbits 

H I in 01 and there is a method to classify them [35]. 

From now on suppose that g is simple. We will say that a 6'-group is exceptional (respec- 
tively classical), if g is exceptional (respectively classical). The classical case allows a more 
or less uniform treatment, since here everything is determined by the canonical embedding 
into an appropriate 0[„, see e.g. [34]. For inner automorphisms of gl„, the nilpotent orbits as 
well as their closures are described by Kempken ([25]). The complete answer, for all classical 
types and all automorphisms, is not known, but there does not seem to be any profound 
difficulty in getting it. 

More interesting representations arise in the context of exceptional 6'-groups. Here several 
orbit classifications were carried out along the lines of [35]. To mention a few, [38], [1], 
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[2], [16]. In these papers, all orbits, not only the nilpotent ones, were described. More 
recently, Pervushin treated one 6'-group in type Ej [28], he also got the closure diagram of 
the nilpotent orbits [29]. 

Despite the possibility to treat each particular exceptional ^-group by hand, the "classical" 
uniformity is lost and one faces a long list of different examples. Dealing with all of them 
by hand is at least difficult. Several computer algorithms for classifying nilpotent Go-orbits 
in 01 have been developed, see [26] and [20]. In this paper, we give a method how to check 
whether a nilpotent orbit Gqx lies in the closure Goy of another nilpotent orbit Goy. 

Each nilpotent element e G gi can be included into an s[2-triple (e, h, f) with h G Qo- Our 
method relies on the fact that this h is also a characteristic of e in the sense of Kempf and 
Hesselink, i.e., it gives rise to a one-dimensional torus in Gq that takes e to zero fastest. 
Another important ingredient is that G^e coincides with a Hesselink stratum, the set of 
all elements in Qi having ^ as a Hesselink characteristic, see [37, Section 5]. Therefore 
Gqc = Go (V>2(/i)), where V>2{h) is the linear span of all vectors f G 0i such that [h, v] = kv 
with k > 2. 

An orbit Gqc' lies in Gqc if and only if its intersection with V>2{h) is non-empty. When 
examining Gqc' fl V>2{h), we replace Gq by the union of its Bruhat cells. Futher, let h' be a 
characteristic of e' and Wq the Weyl group of Gq. Then Propositions 3.1 assures that Gqc' is 
contained in G^e if and only if there is w E Wq such that U{w) = V2{h') fl V>2{'wh) contains 
a point of G^e' (here V2{h') is the set of all vectors v E Qi such that [h', v] = 2v). If this is 
indeed the case, then U{w) fl Goe' is an open dense subset of U{w) and by taking a random 
u G U{w) we can find an element of Gqc' with probability almost one. In order to prove that 
the intersection in question is empty, we compute the dimension of a maximal Z(/i')-orbit 
intersecting U{w) for the centraliser Z{h') C Go of h'. Recall that dim Z{h')v < dim Z{h')e' 
for all elements v in V2{h') \Goe' (see Lemma 2.5). To loop over an orbits of the Weyl group, 
we use its parametrisation as a tree with edges given by simple refiections (Section 4). 
Other tools are described in Sections 2.1, 5, and 6. In particular, to prove a non-inclusion 
Gqc' (^ GqC for some orbits, we use Theorem 2.9, which is a general statement on Z-graded 
reductive Lie algebras and is interesting in itself. It already appeared in the literature and 
was proved by Kac in a particular case [23], see Remark 2.10 for a detailed discussion. 

First examples of ^-groups are provided by the simple Lie algebras themselves, i.e., in the 
case where the automorphism is the identity. Then one asks for the Hasse (closure) diagram 
of the nilpotent orbits in g. The two most difficult, largest exceptional Lie algebras, of types 
Ej and Eg, were treated by Mizuno ([27]). Later his results were verified and corrected by 
Beynon and Spaltenstein ([4]). The implementation of our method in GAP also works for q. 
We have computed the Hasse diagrams for the Lie algebras of exceptional type, and obtained 
the same diagrams as in Spaltenstein's book [32]. 

The same problem for real exceptional Lie algebras has been studied by Djokovic in a 
series of papers [8, 9, 10, 11, 12, 13, 14, 15]. If 0k is a non-compact real form of g and IB C Sr 
is the Lie algebra of a maximal compact subgroup in Gk, then the complexification t(C) 
of Ms a symmetric subalgebra, i.e., 1^(C) = q^ for 6 of order two. The Kostant-Sekiguchi 
correspondence (see e.g. [6, §9.5]) establishes a bijection between nilpotent Gr orbits in 
giR and nilpotent Go-orbits in gi. Moreover, according to [3], this bijection preserves the 
closure ordering. For each automorphism of order 2 of each exceptional complex Lie algebra, 
Djokovic gives the closure diagram for the nilpotent orbits. With the implementation of our 
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method in GAP we have also computed these diagrams. The results of our computations 
were the same as those of Djokovic, except in one case in type Eg. The difference is described 
in Section 8.1. 

The finite order automorphisms of g have been classified by Kac ([24]), up to conjugacy. A 
conjugacy class of automorphisms is identified by its Kac diagram. Here we briefly indicate 
how this works for inner automorphisms, for more information we refer to [36, Chapter 3, 
§3] and [22, Chapter X]. Let $ be the root system of q with a basis {ai, . . . ,ai}. Let ao 
denote the lowest root of $. The Dynkin diagram of the roots ao,ai, . . . ,ai is the extended 
Dynkin diagram of $ (or of g). Let n^ G N be such that ao = — ^j^o*^**^* ^^^ ^^^ '^0 = 1- 
Take / + 1 non-negative integers Sq, . . . ,si with gcd(so, ■ ■ ■ ,si) = 1 and set m = X]j=o'^«'^«- 
Let a; G C be a primitive ?7i-th root of unity. Then a linear map ^ : g — > g that multiplies 
vectors in the root space g^. (0 < i < /) by u'^* uniquely defines an automorphism of g of 
order m. The Kac diagram of this automorphism (or, more precisely, of its conjugacy class) 
is the extended Dynkin diagram with labels sq, . . . ,si. The automorphisms that will appear 
in the examples in this paper all have the labels Si equal to or 1. We will give the Kac 
diagram of such an automorphism by colouring the nodes of the extended Dynkin diagram: 
a black node means that the corresponding label is 1, otherwise it is 0. 

There is also an easy way to read the ^-representation from the Kac diagram of an inner 
6. The group Gq contains a maximal torus of G and the semisimple part of go is generated 
by all root spaces Qai {0 < i < I) with Sj = 0. The lowest weights of Qi (with respect to Go) 
are in one-to-one correspondence with the roots labeled with 1. 

There are two instances of 6 groups, one in Ej and one in Eg, where Go-orbits correspond to 
isomorphisms classes of two-step nilpotent (or metabelian) Lie algebras n such that n' = [n, n] 
is the centre of n and either dim(n/n') < 6,dimn' < 3; or dim(n/n') < 5,dimn' < 5, see 
Section 8.3 and [16]. The nilpotent orbits correspond to those Lie algebras, whose structure 
tensor can be contracted to zero by a unimodular change of coordinates. Here taking closure 
of a nilpotent orbit can be interpreted as the degeneration of the encoded Lie algebra. The 
Lie algebra structures on a given vector space form an afiine algebraic variety and some of 
its properties depend on the degenerations, see e.g. [33]. In the Appendix we present the 
Hasse diagrams for the nilpotent orbits of both these ^-representations. 

We have also computed the closures of the nilpotent orbits of SL9(C) in A'^(C^), see 
Figures 1, 2. This is a ^-representation treated in [38]. 

Section 7 contains a few further observations on algebraic actions. We briefly discuss 
difficulties arising in developing a practical algorithm for describing the closure (Section 7.3); 
outline possible modifications in our algorithm; and present a parametrisation for a set of 
the double cosets of a Weyl group (Section 7.2), which appeared as a byproduct of our 
constructions. 



2. Preliminaries 

In this section we present some results, mainly taken from [37], on which our method is 
based. 

Throughout we let [)o be a fixed Cartan subalgebra of go. The Weyl group of the root 
system of go with respect to i)o will be denoted Wq. We have Wq = Ncoi^o) / Zcoi^o) ■ Hence 
every w G Wq can be lifted to a g E Go such that g\^g = w. Usually we will denote these two 
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elements by the same symbol. The group G is assumed to be simple unless explicitly stated 
to the contrary. 

We say that an s[2-triple (/i, e, /) is homogeneous if e G 0i, /i G go? / ^ 0-i- 

Let us recall a few useful facts (see [34], [35]): 

(1) For a nilpotent element e G gi there exist h E Qq, f E 0_i such that {h,e,f) is an 
s[2-triple. The element h is called a {Dynkin) characteristic of e. 

(2) Let {h',e',f'), {h,e,f) be two homogeneous s[2-triples. Then e' and e are Go- 
conjugate if and only if {h',e',f'), {h,e,f) are Go-conjugate, if and only if h' and 
h are Go-conjugate. 

Thereby a nilpotent orbit Goe corresponds to a unique Go-conjugacy class of homogeneous 
s[2-triples {h,e,f). Also, we may assume that h lies in [)o. Furthermore, after possibly 
replacing hhj a. VFo-conjugate, we may assume that h lies in a fixed Weyl chamber Go of f)o- 
Then h is uniquely determined by the orbit GoC. 

Throughout we will write V for the space gi. Then for h E i)o we set 

Vk{h) = {veV\[h,v] = kv}, V>k{h) = ^Vi{h). 

l>k 

Also we consider the parabolic subalgebra p{h) C Qq, which is the sum of the eigenspaces 
of h with non-negative eigenvalues. Let P{h) denote the connected subgroup of Go with 
Lie algebra p{h). We let ^{h) be the centraliser of h in go- Let l{h) denote the orthogonal 
complement of h in ^{h), with respect to the Killing form of g. Let Z{h) and Z{h) be 
connected subgroups of Go with Lie algebras ^{h) and 3(/i), respectively. 
Now we will borrow two theorems from [37, Section 5]. 

Theorem 2.1 ([37, Theorem 5.4.]). Let e E Qi be nilpotent and nonzero. Let h e\)q he such 
that e G V>2(/i). Then h is a characteristic of e if and only if the projection of e on V2{h) is 
not a nilpotent element with respect to the action of the group Z{h). 

Remark 2.2. In [37, Section 5], the term "characteristic" is used in a different sense, it is 
not necessarily a Dynkin characteristic. However, following the lines of Example 3 in [37, 
Section 5.5], one can show that the orbit Z{h)e is closed in V, if /i is a Dynkin characteristic of 
e. Therefore a Dynking characteristic of e is also a characteristic in the sense of Theorem 2.1. 

The next theorem is the second part of [37, Theorem 5.6.] and Corollary 2.4 is an imme- 
diate consequence of Theorem 2.3. 

Theorem 2.3. Let O = Gqc be a nilpotent orbit in V ; and let h E Qq be a characteristic of 
e. ThenO = Go{V>2{h)). 

Corollary 2.4. Let O' = Gqc' , O = GqC be two nilpotent orbits in V . Let h' , h be Dynkin 
characteristics of e' , e, respectively. Then O' G O if and only if V>2{h) contains a point of 
O'. 

We use the notation g^^^. for the intersection of Qi and the centraliser g^,. C g of x G g. 
Next we have two lemmas that we will use in the sequel. The first one is an immediate 
consequence of Theorem 2.3. 

Lemma 2.5. Let {h,e,f) be a homogeneous sl2-triple. Then Z{h)e is dense in V2{h). 
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Lemma 2.6. Let {h,e, f) be a homogeneous 5i2-triple, and O = GqC. Then h is a (Dynkin) 
characteristic of all elements of O (1 V2{h). 

Proof. In fact, we are going to prove that Y = O D V2{h) is a single Z(h)-OYhit, i.e., this 
intersection is equal to Z{h)e. Let y be an element of Y . Let fj be the eigenspace of h 
in go with eigenvalue i and let fj^y = fi H go,?/ be the centrahser of y in fj. Since Z{h)e is 
dense in V2(/i), the element y lies in its closure. In particular, taking the limits one sees 
that dimfj^j^ > dimfj^e for all i. On the other hand dimgo,e = dimgo,j/, since these are the 
elements of the same Go-orbit. Taking into account that go,j; = @U,y ^^"^ ^^at fo,j; = i{h)y, 
we conclude that the Z(/i)-orbits of e and y have the same dimension. Since Z{h)e is the 
unique Z(/i)-orbit of the maximal dimension, y G Z{h)e. 

Now the statement about characteristics is obvious. D 

2.1. Reduced ^-groups. We conclude Section 2 with a few statements concerning 6'-groups 
appearing from Z-gradings. In this part of the paper, G is an arbitrary (not necessarily 
simple) reductive group. A Z-grading of g = LieG is defined by a diagonalisable one- 
parameter subgroup of Aut(g) and therefore by the eigenvalues of some h E g, i.e., g^ = {■C ^ 
g I [h, C^] = s^}, see e.g. [36, Ch. 3, Sec. 3.3]. Without loss of generality we may assume that 
h G [q,q\- Here all elements of gi are nilpotent and therefore there is a dense open Go-orbit 
in gi. 

Let p : G — 7- GL(iy) be a faithful linear representation of G on a finite- dimensional vector 
space W. We use the same letter p for its differential p : q ^ qK^) ^^^ define a non- 
degenerate G-invariant symmetric scalar product ( , ) on g by setting {x,y) := tr(p(x)p(y)) 
for x,y E Q. Note that the restriction of ( , ) to each non-abelian simple factor of g is the 
Killing form multiplied by a positive rational number. One of the benefits of this choice 
is that {h, h) > 0, whenever h j^ 0, and this is assumed to be the case. More generally, 
(s, s) > for all non-zero s G [g, g] that have rational eigenvalues on g. 

Let go C go be the orthogonal complement of h with respect to ( , ) and Go C Go a 
connected algebraic group with Lie Go = go- Then the action of Go on gi is said to be a 
reduced ^-representation and Go a reduced ^-group. Note that Go = Go(exp(C/i)). 

Lemma 2.7. Let x G gi. Then Gqx = Gqx if and only if [go,x] = [go, a;], and the equality 
takes place if and only if the orbit GqX is conical. 

Proof. If Gqx = Gqx, then clearly [go, x] = [go, x]. Other way around, the equality of tangent 
spaces implies that dimGox = dimGox. Since Go is a normal subgroup of Go, the same 
holds for all elements in Gqx and the two orbits coincide. Finally, being conical means that 
Cx C Gox, or, equivalently, Cx C [go,x]. Therefore Gqx is a conical orbit if and only if 
there is the equality of orbits or their tangent spaces. D 

Lemma 2.8. Suppose that 2h is a Dynkin characteristic of x G gi. Then [go,a;] = Qi, but 
[So,x] T^gi. 

Proof. If [go,x] 7^ gi, then there is v G g_i such that ([go,x],f) = and also (go, [x, u]), 
where [x,v] G go. Since the scalar product is non-degenerate on go, we obtain [x,v] = 0, 
which contradicts the 5l2-theory. 

There is an element y G g_i such that y, 2h, and x form an s[2-triple. For this y we have 
{y, [go,x]) = 0, because (2/i, go) = 0. The inequality follows. D 
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Theorem 2.9. Let G be an arbitrary reductive group and the objects Qi, Gq, Gq as above. 
Suppose that x G 0i. Then GqX 7^ GqX if and only if 2h is a Dynkin characteristic of . 



X. 



Proof. Let h & Qq he a. Dynkin characteristic of x. We can write it a.s h = ah + h^ with 
a G C and h^ G go- Since [/i, x] = 2x = ax + [hQ.,x\, either Cx C [floj^;] or a = 2 and 
[/lo,^:] = 0. In the latter case ih^ho) = 0. Taking into account the equahty {h,hQ) = 0, we 
get that {ho, ho) = 0. Since /i is a Dynkin characteristic, it lies in [q,q]. Hence /iq G [0,g], 
because h also does. Moreover, eigenvalues of ad(/i) are integers by the construction, and the 
same holds for ad(/i), because it conies from an s[2-triple. Since [h,h] = 0, the eigenvalues 
of ad(/io) are integers as well. According to our choice of the scalar product, the equality 
{ho, ho) = is possible only if /iq = 0. One concludes that 2h is a Dynkin characteristic of x. 

We have shown that if 2h is not a Dynkin characteristic of x, then [ho, x] = bx with 6 G C , 
in particular, Gox is a conical orbit. By Lemma 2.7, Gqx = Gox. 

If 2h is a Dynkin characteristic of x, then [qo,x] 7^ [Qo,x] by Lemma 2.8 and therefore 

GqX 7^ GqX. D 

Remark 2.10. In case q is simple and the representation of Go on gi is irreducible. The- 
orem 2.9 was proved by V. Kac, see [23, Proposition 3.2]. It is also mentioned without a 
proof in [37, Section 8.5] that the statement holds for an arbitrary reduced ^-group. Since 
we could not find a general case proof in the literature, we decided to include it here. 

3. Criteria for inclusion 

In this section we state and prove the main criterion (Proposition 3.1) that we use for 
deciding whether a given nilpotent orbit is contained in the closure of another given nilpotent 
orbit. This reduces the problem of checking inclusion to a finite number of checks, each 
corresponding to an element of a certain orbit of the Weyl group Wo- Subsequently we give 
some observations that help when using the criterion. 

Proposition 3.1. Let the notation be as in Corollary 2.4. Then O' G O if and only if there 
is a w E Wo such that U = V2{h') fl V>2{wh) contains a point of O' . Moreover, in that case 
the intersection of U and O' is dense in U. 

Proof. The "if" part follows directly from Theorem 2.3. Therefore suppose that O' C O. 
By the Bruhat decomposition we have that 

Go= U P{h')wP{h). 

w&Wo 

By Theorem 2.3, 

= Go{y>2{h)) 

= IJ P{h')wP{h){V>2{h)) 

= U P{h')w{V>Ah)). 

Let {h' , e! , /') be a homogeneous s[2-triple. Then it follows from the above that there exist 
w G Wo, p G P{h'), and x G V>2{h) with e' = pwx, or, equivalently, p~^e' = wx. 



CLOSURES OF NILPOTENT ORBITS 7 

Next P{h') = Z{h') x A^, where A^ is a connected subgroup of Gq whose Lie algebra is the 
sum of the eigenspaces of h' in go with positive eigenvalues. In particular, p~^ = In with 
/ G Z{h') and n & N. Since e' G V2{h'), we have ne' = e' + y, where y G V>3(/i'). Now 
p~^e' = le' + ly, with le' G V2{h') and ly G V>3(/i'). In particular, p~^e' lies in V>2(/i')- Since 
p~^e' = wx and wx G V>2(tf/i), it also lies in U = V>2(^') H y>2(tf/i). 

The elements h' and w/i commute and thereby U is stable under the action of h'. That is, 
U is the direct sum of /I'-eigenspaces. It follows that U contains le', which is obviously an 
element of O'. Moreover, le' G V2{h') and hence le' G U, where U = V2{h') n V>2{wh). 

By Theorem 2.1, an element v & U lies in O' if and only if it is not nilpotent with respect 
to the action of Z{h'). Threfore if the intersection of U and O' is not empty, then it has to 
be open and dense. D 

Proposition 3.2. Let {h',e',f'), {h,e,f) be homogeneous sl2-triples, with e' G O' , e G (9. 
Let K denote the Killing form of g. If K,{h',h) < K,{h',h') then V2{h') fl V>2(/i) contains no 
points of O' . 

Proof. Note that h G ^{h'), hence h = ah' + 1, where a G C and t G Jih'). Moreover, 

_ n{h',h) 

"~ K{h',h')' 

which is in Q and < 1. Hence t has only positive eigenvalues on V2{h') fl V>2(/i). Let T 
be the connected subgroup of Gq whose Lie algebra is spanned by t. Then all elements of 
V2{h') n V>2{h) are nilpotent with respect to T, and in particular with respect to Z{h'). 
Hence by Theorem 2.1 and Lemma 2.6, the former space contains no points of O' . D 

Let [ be a Lie algebra acting on a vector space M. Then for t; G M we denote its stabiliser 
by [,;, i.e., 

ti, = {x G I I x-v = 0}. 
The set of f G M with dim ly minimal is open and dense in M. 

Proposition 3.3. Let {h',e',f'), {h,e,f) be homogeneous sl2-triples, with e' G O' , e E O. 
Let d = diin^{h')e>- Let d' be the minimal dimension of ^{h')u, for u G V2{h') fl V>2(/i). Then 
d < d' . Moreover, V2{h') fl V>2{h) contains a point of O' if and only if d = d' 

Proof. By Lemma 2.5, the stabiliser of e' in ^{h') has minimal possibile dimension. 

Furthermore, ii d = d' then there is u G V2{h') fl V>2{h) such that dim3(/i')u = dim.^{h')e'. 
Hence the dimension of the Z{h')-oihit of u is the same as the dimension of Z{h')e' . So 
Z{h')u is dense in V2{h') as well. The conclusion is that Z{h')e' = Z{h')u, and u lies in 
O'. D 

Proposition 3.4. Let {h',e',f'), {h,e,f) be homogeneous 5l2-triples, with e E O, e' E O' . 
Set U = V2{h') n V>2{h). Let n = Ng^{U) = {x G flo I [x, U] C U}. Let u G U ; if [n, u] = U, 
and u ^ O' , then U has no point of O' . 

Proof. Indeed, if U has a point of O' , then the intersection of O' and U is dense in U . 
But also the A'^G'y(f/)-orbit of u is dense in U . So the two sets must intersect, which is not 
possible. D 
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4. Orbits of the Weyl group 

In our algorithm we need to loop over an orbit Woh, where h G f)o- In this section we 
briefly describe how this is done. For simplicity we assume that the centre of Qq is zero. If 
this is not the case then Qq has to be replaced by its derived subalgebra [007 00]; and f)o by 
its intersection with [flcflo]- 

We let K denote the Killing form of g. Since it is non-degenerate on i)o it gives an isomor- 
phism 1)q — )■ f)o, a !-!■ a. This yields an inner product on [)q by (a, (3) = ^{a, f3). 

Let $0 be the root system of go with respect to i)o. Let Aq = {ai, . . . , a/} be a fixed basis 
of <l>o. The corresponding set of positive roots will be denoted ^q. 

For a G $0 we set 

V 2d 

a = G ()o- 

[a, a) 

The Weyl group I^o is generated by the simple reflections Sj = s^^. For h ^ i)^ we have 

Si{h) = h — ai{h)a^ . 

We use a basis hi, . . . ,hi of l)o, defined by ai{hj) = 6ij. Then, ii h = J2i (^ihi, we get 
Sj{h) = h — ajOj . The elements h of which we compute the IVo-orbit, lie in an s[2-triple. 
This implies that the coefficients of h with respect to this basis are integers. The dominant 
Weyl chamber Cq consists of the elements of f)o having non-negative coefficients with respect 
to the basis hi, . . . ,hi. 

Now let /i G f)o be the element of which we want to compute the orbit Woh. Since every 
orbit of Wq has a unique point in Co, we may assume that h E Cq. Let h G W^h, then we 
define the length of h, denoted l{h), as the length of a shortest w G Wq with h = wh. Then 

£(h) = |{aG<l>^ I aCh) < 0}|. 

This implies that i{sih) = i{h) + 1 if and only if a^ > 0, where h = ^^aihi. We use a 
criterion due to Snow ([31]): 

Lemma 4.1. Let h = ^^ Oj/ij be an element ofWoh of length k + 1. Then there is a unique 
h of length k in Woh such that 

• there is a simple reflection Si with Si{h) = h, 

• ttj > for i < j < L 

Let h, h be as in the previous lemma. Then we say that h is the predecessor of h, and 
conversely, that h is a. successor of h. Let h = Y2- hihi be a given element of W^h of length 
k. Then it is straightforward to determine its successors. Indeed, let i be such that hi > 0, 
and write Si{h) = ^^ Oihi. Then this element is of length k + 1, and it is a successor of h if 
and only if aj > ioi i < j < I. 

This means that we can define a tree: the nodes are the elements of Woh, and there is 
an edge from h to h ii and only ii h is a. successor of h. By traversing this tree, we can 
efficiently loop over Woh. Every element of Woh comes at the cost of applying one reflection. 
Moreover, we do not obtain the same element of Woh twice. 

Remark 4.2. We finish this section with an observation that will be used later. Let h' be 
an element of Co- Let h G Woh be of length k and suppose that Sjh is of length k + 1. Write 
h = ^^aihi; then, as seen above, Oj > 0. Hence 



K 



{sjh, h') = K{h, h') — ajhi{a^j , h') < k(/i, h'). 
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Furthermore, equality happens if and only if aj{h') = 0, which is equivalent to Sj lying in 
the stabiliser of h'. 

5. Complement of the dense orbit 

According to Proposition 3.1, we will have to check whether a subspace U C V2{h) contains 
a point of the dense orbit Z{h)e. If U contains a point of Z{h)e then the intersection of U 
and Z{h)e is dense in U (Proposition 3.1). So in that case, by trying random elements of 
t/, we quickly find a. u &U lying in Z{h)e; thus proving that the intersection is non-empty. 
The most difficult part of the problem is to prove that U contains no points of Z{h)e. Here 
we present two possible solutions. 

Let Vi, . . . ,Vs and Xi, . . . , x„ be bases of V2(h) and ^(h) respectively. Let also Wi, . . . ,Ws 
(with Wi G V2{h)*) be the dual basis. Let B denote the action matrix for the representation 
of 3(/i) on V2{h)*. To be more explicit, the entries of B are elements of V2{h)*, bij = Xi-Wj. 
For V & V, let By denote the restriction of -B to f . The entries of this new matrix are 
[xi-Wj]{v) = Wj{[v,Xi]). In the same spirit, we can define the restriction of B to U, Bu, to 
be a matrix with entries in U*. The rank of Bu is calculated over the field C(f/) (note that 
U* C C([/)). 

Using the fact that [C,, v] = (with C, G ^{h)) if and only if Wi{[C,, v]) = for all i, one can 
easily deduce that 

(i) dim.^{h)y = n — ranki?^ for all v G V2{h)] 
(ii) dim Z{h)v = ranki?^; 
V / (iii) maxdimZ(/i)n = ranki?{7; 

uGU 

(iv) U n Z{h)e 7^ if and only if iai\kBu = s. 

Depending on s and n, computing the rank of Bu over a function field may turn out to 
be rather time consuming. For this reason we also consider an alternative method, based on 
another characterisation of the elements in V2{h) \ Z{h)e, which comes from Theorem 2.9. 

Proposition 5.1. Take v G V2{h). Then the three conditions: Z{h)v = Z{h)v, [i{h)^v\ = 
[1(h), v], and v G V2{h) \ Z{h)e, are equivalent. 

Proof. We are going to identify Z{h) with a reduced 6'-group. To this end, for each i G Z, 
set i := i mod m, if 6 has a finite order m; and i := i otherwise. Then [ = ^ Ij, where 

h = {9i)2i{h), is a Z-graded Lie subalgebra of g with [i = V2{h) and Iq = 3(^)- Let L C G 
be a connected subgroup with Lie L = I. Since k defines a non-degenerate pairing between 
{Qi)2i{h) and (gj)_2i(/i) with j = —i, we get a non-degenerate L-invariant scalar product 
( , ) := k|[ on [. In particular, [ is a reductive subalgebra. Here e G (0i)2(^) = ti, / G 
(g_i)_2(/i) = Li and therefore h G [I, I]. Note that the Z-grading on [ is defined by the 
eigenvalues of h/2. 

Recall that q is assumed to be simple. Restricting the adgoint action of G to L we get a 
faithful representation p of L on g such that (x, y) = tr {p{x)p{y)) for x,y E I and { , ) = k\[. 
Thus, we are in the setting of Section 2.1 and can apply Theorem 2.9 to the Z-graded 
reductive Lie algebra [. Here Lq = Z{h) and li = (gi)2(^) = V2{h). 

We have [l{h),v\ = [l{h),v] if and only if the Z(/i)-orbit Z{h)v is conical. Besides, h is 
a Dynkin characteristic of all elements in Z{h)e. Therefore both equivalences follow from 
Theorem 2.9. D 
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Assume that the basis Xi, . . . , x„ is chosen in such a way that Xi = h and X2, . . . , x„ form 
a basis of l{h). Let B he a. submatrix of B consisting of the last n — 1 rows (corresponding to 
the Lie subalgebra ]){h)). Let also Bjj be the restriction of B to U. Since dim Z{h)v = ranki?„ 
and dim Z{h)v = ranki?^, Proposition 5.1 gives us the following: 

(5.2) U n Z{h)e = if and only if iai\kBu = ranki?f/. 

In other words, either ranki?^ = lankBu or u is an element of Z{h)e and ranki?„ = s. The 
equality in equation (5.2) is satisfied if and only if the first row of Bu lies in the linear span 
of the rows of Bu. In order to check this we use the following steps. 

(1) Take a random u E U, compute the rank of Bu, say ranki?^ = r. 

(2) Find an r x r non-zero minor of Bu, without loss of generality suppose that it is given 
by the first r rows and the first r columns. 

(3) Check whether the first row of Bjj is contained in the span of the first r rows of Bir. 

If u G f/ is generic, i.e., ranki?^ = TankBu, then the first r rows of Bu span the row space 
of Bu- Hence step (3) verifies whether the first row of Bu is contained in the row space of 
Bu- Moreover, this will be the case if and only if f/ fl Z(/i)e is empty. (Also note that the 
check in the third step can be done by computing s—r minors of size r+1.) Even if u is not 
a generic element, it may still be true that the first row of Bu is contained in the span of 
the first r rows of Bu, and the above procedure will prove that U fl Z{h)e = 0. 

In many cases it is easier to carry out this procedure than to check the inequality rankBu < 
s. For example, some 32 x 38-matrices Bu appeared while checking non-inclusions for a half- 
spin representation of Dg (line 3 in Table 1) and 2760681 minors would have to be computed 
for them. In other cases it may be easier to deal with the whole matrix, if, for example, Bu 
contains a zero column. 

It is not obvious beforehand which choice is the best. In the implementation of our 
algorithm we do the following: ii s — n < s — r, then it is checked whether lankBu < s. 
Otherwise we check whether the first row if Bu is contained in the first r rows of Bu, using 
the procedure outlined above. We do not claim that this always gives the best choice, but 
some choice is better than none. 

If it turns out that the first row of Bu is not contained in the span of the first r rows 
of Bu, then it may still be the case that the intersection is empty (if this happens, then 
necessarily lankBu > ranki?„). Then we will have to compute the rank of Bu- However, 
the probability of this event can be made arbitrarily small. 

6. The main algorithm 

Here we describe our algorithm for deciding whether one of the two given nilpotent Go- 
orbits in Qi lies in the closure of the other. 

First we consider the following problem: given a homogeneous s[2-triple {h, e, /) and 
e' G ^(/i), decide whether e' G GqC- We have a straightforward solution for that, based on 
Lemma 2.6. The existence of /' G g_i with [h, /'] = —2/' and [e', /'] = h is equivalent to 
a system of linear equations. We solve this system; if it has a solution then e' lies in GqC, 
otherwise it does not. 

Throughout we fix a basis of the root system of Qq with respect to P)o. Then the Weyl group 
Wq is generated by the reflections corresponding to the elements of this basis. Furthermore, 
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this choice also fixes a dominant Weyl chamber Cq G i)o. As before we let k, denote the 
Killing form of g. 

Algorithm 1. Input: two homogeneous 5l2-triples, {h', e', /'), {h, e, /), such that h', h G Cq. 
Output: TRUE if O' = Gqc' is contained in the closure oi O = Gqc, false otherwise. 

(1) If K{h', h) < K,{h', h') then return FALSE. Else go to the next step. 

(2) For all elements wh G Woh do the following: 

(a) If K{h',wh) > K{h\h') then: 

(i) Select a random u eU = V2{h') fl V>2(w/i). 
(ii) If M G O' then return TRUE. Otherwise go to the next step, 
(iii) Set n = Ng^^{U). If [n,u] ^ U then decide whether U fl Z{h)e is empty 
using the methods of Section 5. If the intersection is not empty then return 

TRUE. 

(3) If in the previous loop true was never returned, then return FALSE. 

Proposition 6.1. The previous algorithm terminates correctly. 

Proof. It is obvious that the algorithm terminates, we must show that the output is correct. 
We claim that the algorithm checks whether there is a wh G W^h such that U{wh) = 
V2{h') n V>2{'wh) contains a point of O'. Then by Proposition 3.1 the output is correct. 

First of all we note that since h',h G Cq we have that the maximal value of K{h',wh), 
for wh G Wh, is K{h',h) (see Remark 4.2). Therefore, if K{h',h) < K{h',h') then no space 
U{wh) contains a point of O' (Proposition 3.2). So in this case we are immediately done. 

Otherwise we inspect every wh G Woh. If K{h', wh) < k(/i', h') then U{wh) contains no 
points of O' by Proposition 3.2. So then we can discard it. Otherwise we select a random 
u G U{wh). If n G O', then we are done. If not, and [n,u] = U{wh), then U{wh) has 
no points of O' by Proposition 3.4. Finally U{wh) contains an element with a minimal 
dimensional stabiliser if and only if U{wh) has a point of O' by Proposition 3.3. D 

Remark 6.2. Note that if U{wh) contains a point of O', then the set of such points is open 
and dense in U{wh), by Proposition 3.1. Hence in that case the random choice has a high 
probability of finding an element u ^ O'. 

Remark 6.3. Now we make some observations that help to execute the algorithm more 
efficiently. 

• Of course, we apply the algorithm only if dim O' < dim O, as otherwise there is no 
inclusion. 

• Inclusion also implies that dimQk^e' > dimg^^e for all k; so if that condition is not 
fulfilled, we also do not apply the algorithm. 

• When looping over the orbit Woh we use the tree structure described in the previous 
section. When doing this, several shortcuts can be made. First of all if K{h', wh) < 
K{h', h') then the entire subtree below wh can be discarded, by Remark 4.2. Second, 
if l^(/i') n V>2{'wh) contains no points of O', and the successor of wh is Siwh, where 
Si lies in the stabiliser of h', then also V2 (/?-') fl V>2{siwh) contains no points of O'. 
Hence in that case we can immediately jump to the next element of the orbit. 

• We collect the subspaces U{wh) = V2{h') nV>2{wh) that appear during the execution 
of the algorithm. If a certain such subspace is contained in one that was treated 
before, then we already know that it contains no points of O'. So in that case we can 
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immediately go to the next round. All calculations are done in a basis Vi, . . . ,Vs of 
V2{h') consisting of f)o-eigenvectors. Each subspace U{wh) is a linear span of Vi such 
that {wh)-Vi = ttiVi with Oj > 2. Therefore storing and verifying inclusions among 
the U{wh) is a binary problem. 

We have implemented this algorithm in the language of the computer algebra system GAP4 
([17]), on top of the SLA package ([19]), which has implementations of algorithms to list the 
nilpotent orbits of a ^-group. One of the main problems for the practical computation lies 
in the methods of Section 5, where minors of a matrix with entries in a function field have 
to be computed. For these computations we use the computer algebra system Magma ([5]). 
We have chosen this system, because it has very efficient implementations of algorithms to 
compute the determinant of a matrix with entries in a function field. 

In Table 1 we collect some experimental data with respect to the implementation of our 
algorithm. All computations have been performed on a 3.16 GHz machine. 



Kac diagram of d 






^ orbits 



69 



36 



115 



101 



105 



GAP 



1003 



124 



900 



444 



335 



Magma 



397 



0.18 



0.09 



0.04 



Table 1. Running times (in seconds) for the algorithm applied to several 
automorphisms of the Lie algebra of type Eg. The ffist column lists the order 
of 6, and the second column its Kac diagram. The third column has the 
number of nilpotent orbits. The fourth and fifth columns display, respectively, 
the time needed for executing the GAP part of the program, and the time spent 
in the Magma part. 



From the table we see that the GAP part essentially has no problems, also with large 
examples. On some occasions it is not necessary to execute the Magma part, as with the 
second automorphism in the table. On other occasions this part has a trivial running time, 
as with the last three examples. However, it also happens that a fair amount of time is spent 
in the Magma part, as with the ffist example. 
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7. Further remarks on groups and orbits 

Here we collect some theoretical observations that could have been used in the algorithm, 
but turned out not to be of much practical value. 

Let (/, h, e) be a homogeneous sl2-triple as before, we also keep all the previous notation, 
including V2{h). First we consider the actions of Z{h) and Z{h) on V2{h) more closely. As 
was already mentioned, Z{h) acts on V2{h) with a dense open orbit, Z{h)e. 

7.1. The semi-invariant P. The stabiliser 3(/i)e, being the centraliser of an s[2-subalgebra 
generated by e, /i, and /, is reductive and therefore the orbit Z{h)e is an afiine space. This 
implies that the complement V2{h) \ Z{h)e is a divisor and is a zero-set of a single semi- 
invariant polynomial, say P. To check whether a subspace U = U{wh') = V2{h) fl V>2{u!h') 
intersects the dense orbit, one just has to look on the restriction of P to U. In this terms, O 
lies in the closure of O' if and only if there is U{wh') such that P is non-zero on it. This could 
be a replacement for both: choosing a random element in U and generic rank considerations. 

One can try to compute a polynomial P by hand. This may involve typing errors and 
time-consuming calculations. It is also possible to get P from the matrix B with entries 
bij = Xi-Wj, where {xj} is a basis of ^{h) and {wj} is a basis of V2{h)* (this matrix was 
already considered in Section 5.) The polynomial P is the greatest common divisor of the 
largest, dimV2(/i)x dimV2(/i), minors of B. In some cases the resulting formula is rather 
bulky and not easy to deal with, in some other MAGMA was unable to finish the calculation. 
It turns out that MAGMA checks much more easily that the restriction Bu of 5 to f/ does 
not have the maximal rank, dimV2(/i), than it computes the greatest common divisor of 
minors. Thus we gave up the idea of using P. 

7.2. Double cosets of Weyl groups. Suppose that we have two characteristics h and h' 
lying in the dominant chamber of Wq. Parametrisation of Woh involves a certain numbering 
of simple roots ai, . . . ,ai, see Section 4. This numbering can be arbitrary. The stabiliser 
Wo^h' is a Weyl subgroup generated by s, with aiih') = 0. Assume that rankVFo,/i' = ^ and 
the simple roots orthogonal to h' have numbers from l—r+1 to /. 

Lemma 7.1. Keep the above notation and enumeration of simple roots. Let T be a tree 
parametrising the orbit W^h, constructed according to the principles of Lemma 4-1- Then 
the nodes h and Sih of T with i < l—r are in one-to-one correspondence with the double 
cosets Wo^h'\WQ/WQ^h. 

Proof. First note that h or Sih with i < l—r lies in the dominant chamber of Wo.ft'- Secondly, 
an element Sih with l—r < i does not bring a new double coset, because here Si lies in 
Wo,h'. □ 

This, of course, is not a very effective way for listing the double cosets as the whole orbit 
Wok has to be constructed. However, if some time consuming calculation has to be performed 
for representatives of double cosets, such a treatment may be useful. 

In our situation, collecting subspaces U{wh) turned out to be much more effective than 
refining the Weyl-group tree. The explanation is that one and the same U = U{wh) arises 
for many different elements wh G Wok. 
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7.3. Other algebro-geometric methods. Our algorithm is designed for ^-groups and 
works quite well. There are known some other, more general, approaches, which unfortu- 
nately have a rather small range of application. 

To begin with, consider a linear action a : Q x V ^ V oi an affine algebraic group Q 
on a vector space V. For x & V, the map g ^^ gx from Q to V is regular, i.e., given by 
polynomials. Let I{Qx) < C[V] be the ideal of Qx, i.e., a set of all polinomials vanishing 
on Qx. In [7], algorithmic methods are described for computing generators of the vanishing 
ideal of the image of a regular map. In particular, this can be applied to I{Qx). Here 
I{Qx) equals (a*)"^(/(Qx{x})), where I{Qx{x}) < C[GL(V)x{x}] is the defining ideal of 
the product of the image of Q in GL(\^) and the point x. Once the generators of I{Qx) are 
known, it is straightforward to decide whether a point (and hence the orbit of that point) 
lies in Qx. 

In order to use the algorithms for getting generators of I{Qx), we need as input the 
polynomials defining Q as a subgroup of GL(V^). In our setting, V = Qi and Q is the image 
in GL(y) of Gq, acting on V. In order to get equations for Q, methods from [18] can be 
used. However, both the algorithm for obtaining the polynomials defining Q and the one 
for computing I{Qx) heavily rely on Grobner basis computations. These are extremely time 
consuming. For this reason this method is only applicable to very small examples (e.g., when 
the semisimple part of go is of type Ai, and gi is of dimension 5). 

Now let G C GL[l^] be reductive. Then in the above considerations G can be replaced 
by its big open cell, BwB, where w E W is the longest element in the Weyl group W of 
G and -B C G is a Borel subgroup. In [30], V.L.Popov suggested an algorithm, based on 
this observation, for deciding whether Gy lies in Gx. That algorithm uses a system of linear 
equations in '^("^_]" ) variables, where n = dim^ and d is the degree of G as a subvariety 
of l-i{V). For an irreducible 5-dimensional representation of SL2(C), the number of variables 
equals 56794400. This makes it difficult, if not impossible, to use the algorithm for practical 
computations. 

8. Examples 

In this section we show the output of our programs on several examples. Here we describe 
the examples; the next section contains the Hasse diagrams that we computed with the 
algorithm, as well as tables giving the characteristics of the nilpotent orbits. 

8.1. Symmetric pairs. As was mentioned in the Introduction, the order two case, or, in 
other words, the symmetric case, was studied by Djokovic, because of its relationship with 
simple real Lie algebras. We have checked all symmetric pairs arising from the exceptional 
Lie algebras. The result is that Djokovic diagrams are basically correct, if one takes into 
account the necessary alteration that he found himself, [13], [15]. Our calculations confirm 
these corrections. Apart from this, there are two inclusions missing for one automorphism in 
type Eg. For the involution in question, go is of type Dg and gi is a half-spin representation, 
the corresponding Kac diagram is the third one in Table 1. 

According to our calculations, in Table 2 of [14], 59 -> 53, and 95 — )■ 92 should be added 
(notation as in the mentioned paper: 59 — >■ 53 means that orbit number 53 is contained in 
the closure of orbit number 59). This then results in several other changes. For example, 
[14, Table 2] states "99 — )■ 92, 94, 95" and 92 has to be removed from here, because the orbit 
number 92 does not give rise to an irreducible component of the boundary dOgg. 
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Neither (59,53) nor (95,92) appears in the hst [14, Table 5] of the critical pairs, i.e., pairs 
(a, b) such that the non-inclusion a -/^ b is (or has to be) proved. In both cases, our program 
immediately found that the space U = V2{h') fl V>2{h) contains a point of O'} 

8.2. Trivectors of a nine dimensional space. In [38], the orbits of SL9(C) acting on 
A^(C^) were obtained. This is known as the classification of the trivectors of a nine dimen- 
sional space. The orbits were obtained by realising this representation as a ^-representation. 
Here 6 is an automorphism of order 3 of the Lie algebra of type Eq, with Kac diagram 




We have that go = slg{C) and gi ^ a3(C9). 

In Table 4 we list the characteristics of the (nonzero) nilpotent orbits. A characteristic 
h is given by the values of ai{h), where {ai, . . . , ag} is a basis of the root system of go. 
Moreover, all the characteristics h lie in the dominant Weyl chamber with respect to this 
basis. 

Figures 1, 2 contain, respectively, the top half and the bottom half of the Hasse diagram. 

8.3. The classification of metabelian Lie algebras. A finite-dimensional Lie algebra 
L is said to be metabelian (or two-step nilpotent), if [L, [L,L]] = 0. In [16], Galitski and 
Timashev described the Go-orbits for the two particular ^-representations in order to obtain 
the classification of the metabelian Lie algebras of dimensions up to 9 (over algebraically 
closed fields of characteristic 0). With our algorithm we have computed the closure diagram 
of the nilpotent orbits in both these cases. 

For the first 6'-group, 6 is an automorphism of order 5 of the Lie algebra of type Eg, with 
Kac diagram 




Characteristics of the nilpotent orbits are given in Table 5. The closure diagram is displayed 
in Figures 3, 4. 

For the second 6'-group, 6 is the automorphism of order 3 of the Lie algebra of type Ej, 
with Kac diagram 



o — o 




Characteristics of the nilpotent orbits are given in Table 6. The closure diagram is displayed 
in Figures 5, 6. 

As was already mentioned, in [16] the orbits of these particular ^-groups are used for a 
classification of metabelian Lie algebras. Every orbit corresponds to one such Lie algebra 
(up to isomorphism). Every metabelian Lie algebra L has a signature, that is a pair {m,n) 
where m = dimL/[L,L] and n = dim[L,L]. Let Z <\ L he a maximal Abelian ideal such 
that Z n [L, L] = 0. Then L = L/Z © Z. In the closure diagrams we indicate the signature 
of L/Z as it was computed in [16]. Mostly this is done by writing the label of the node in 
a particular font, according to Tables 2 (for Figures 3, 4) and 3 (for Figures 5, 6). For the 



This can be even verified by hand, would someone wish to do so. 
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orbits corresponding to signatures not present in these tables we have put the signature in 
the diagram, next to the node. 



font 


example 


signature 


ronian 


10 


(5,5) 


bold face 


14 


(5,4) 


italics 


70 


(5,3) 


underhne 


82 


(4,4) 


typewriter 


92 


(4,3) 



Table 2. Fonts for the Hasse diagram in Figures 3, 4. 



font 


example 


signature 


roman 


10 


(6,3) 


bold face 


35 


(6,2) 


italics 


41 


(5,3) 


typewriter 


68 


(4,3) 


overline 


64 


(4,2) 



Table 3. Fonts for the Hasse diagram in Figures 5, 6. 



Taking the closure of a given nilpotent orbit corresponds to the degeneration of the encoded 
two-step nilpotent Lie algebra. Let us explain this. Let VT be a vector space (over C). Then 
a Lie bracket on W can be seen as an element of Hom(A^VF, VT). The group GL(iy) acts 
on Hom(A^iy, ly), and the orbits of this action are in one-to-one correspondence with the 
isomorphism classes of Lie algebra structures on W. Let X, fi E B.oin{A'^W,W); if fi is 
contained in the closure of the GL(H^)-orbit of A, then fj, is said to be a degeneration of A. 
We refer to [21] for an introduction into this concept. 

In relation to the variety of metabelian Lie algebras one considers two vector spaces U and 
V. A metabelian Lie bracket on [/ © ^ is viewed as an element of Hom(A^t/, V). Let L be 
the Lie algebra defined by such an element; then [L, L] C V. The group GL{U)xGL{V) acts 
on Hom(A^f/, V). Two metabelian Lie algebra structures on U (B V are isomorphic if and 
only if the corresponding elements of Hom(A^[/, V) lie in the same GL(f/)xGL(V^)-orbit. 

Write m = dim U, n = dim V, and let W he a. vector space of dimension m + n. If 
A G Hom(A^VF, W) is a metabelian Lie bracket, defining a Lie algebra L on W^ of a sig- 
nature {m,n), then by setting V = [L,L] and taking t/ to be a complement of V in 
W, we get an element of Hom(A^f/, V). This construction preserves isomorphism, i.e., a 
GL(14^)-orbit of metabelian Lie brackets in B.oni{A^U,V) is mapped to a GL{U)xGL{V)- 
orbit in Hom(A^?7, V). Thus metabelian Lie algebras of signature (m, n) are classified by the 
GL(t/)xGL(y)-orbits in B.om.{A^U,V). Moreover, degenerations of these Lie algebras are 
given by the orbits closures. 

Two instances of 6'-groups considered in [16] correspond to the signatures (5, 5) and (6, 3). 
Strictly speaking, the group Go is semisimple in both cases. Therefore, if ^ Gqx with x E Qi, 
then the one-parameter family of Go-orbits Go{ax) with a e C gives only one isomorphism 
class of metabelian Lie algebras. For the nilpotent Go-orbits, there is no difference between 
Goxand {GL{U)xGL{V))x. 
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Note that the quotients LjZ with signatures being smaller than or equal to (5, 3) appear 
in both ^-groups. Therefore, the lower parts of both Hasse diagrams (Figures 3, 4; 5, 6) are 
the same. 

The afiine variety of metabelian Lie algebras structures with signature (5, 5) on a ten- 
dimensional vector space W is irreducible, because of the equivalence with 6'-group orbits. 
According to [16], it has a one-parameter family of the maximal GL(H^)-orbits. The same 
holds for the second case, where signatures are (6, 3) and dim ly = 9, only there is a two- 
parameter family of the maximal orbits [16]. Since isomorphism classes of metabelian Lie 
algebras with smaller signatures are parametrised by nilpotent orbits, there are only finitely 
many of them. In addition, our diagrams show that for each signature (m, ra), where either 
m < 5, n < 5 and {m,n) ^ (5,5), or m < 6, n < 3 and {m,n) ^ (6,3), the afiine variety 
of metabelian Lie algebra structures with signature (m^ n) on W with dim W = m ^ n y& 
irreducible. For example, all metabelian Lie algebras with signature (5, 4) are degenerations 
of LjZ (including LjZ itself), where L corresponds to the orbit number 14 of the first ^-group 
(in E'g); up to isomorphism, see Figure 3. The metabelian Lie algebras with signature (5,3) 
are degenerations of LjZ^ where L is encoded by the orbit number 62 of the first 6'-group as 
well as by the orbit number 27 of the second 6'-group, see Figures 3, 5. 

References 

[1] L.V. Antonyan, Classification of four- vectors of an cight-dimcnsional space, Trudy S em. Vektor. Tenzor. 

Anal, 20(1981), 144-161. 
[2] L.V. Antonyan and A.G. Elashvili, Classification of spinors of dimension sixteen, Trudy Thiliss. Mat. 

Inst. Razmadze Akad. Nauk Gruzin. SSR, 70(1982), 5-23. 
[3] D. Barbasch and M.R. Scpanski. Closure Ordering and the Kostant-Sekiguchi Correspondence. Proc. 

Amer. Math. Soc, 126(1):311-317, 1998. 
[4] W.M. Beynon and N. Spaltenstein. Green functions of finite Chevalley groups of type En{n = 6, 7,8). 

J. Algebra, 88(1984), 584-614. 
[5] W. Bosma, J. Cannon, and C. Playoust. The Magma algebra system. I. The user language. J. Symbolic 

Comput., 24(3-4):235-265, 1997. Computational algebra and number theory (London, 1993). 
[6] D. H. Collingwood and W. M. McGovern. Nilpotent orbits in semisimple Lie algebras. Van Nostrand 

Reinhold Mathematics Series. Van Nostrand Reinhold Co., New York, 1993. 
[7] D. Cox, J. Little, and D. O'Shea. Ideals, Varieties and Algorithms: An Introduction to Computational 

Algebraic Geometry and Gommutative Algebra. Springer Verlag, New York, Heidelberg, Berlin, 1992. 
[8] Dragomir Z. Dokovic. The closure diagrams for nilpotent orbits of real forms of F/^ and G2. J. Lie 

Theory, 10(2):491-510, 2000. 
[9] Dragomir Z. Dokovic. The closure diagram for nilpotent orbits of the real form EIX of E%. Asian J. 

Math., 5(3):561-584, 2001. 
[10] Dragomir Z. Dokovic. The closure diagram for nilpotent orbits of the split real form of E-;. Represent. 

Theory, 5:284-316 (electronic), 2001. 
[11] Dragomir Z. Dokovic. The closure diagrams for nilpotent orbits of real forms of _Bg. J. Lie Theory, 

11(2):381-413, 2001. 
[12] Dragomir Z. Dokovic. The closure diagrams for nilpotent orbits of the real forms E VI and E VII of Ej. 

Represent. Theory, 5:17-42 (electronic), 2001. 
[13] Dragomir Z. Dokovic. Correction to: "The closure diagrams for nilpotent orbits of the real forms E VI 

and E VII of £^7" [Represent. Theory 5 (2001), 17-42.]. Represent. Theory, 5:503 (electronic), 2001. 
[14] Dragomir Z. Dokovic. The closure diagram for nilpotent orbits of the split real form of E^. Gent. Eur. 

J. Math., l(4):573-643 (electronic), 2003. 
[15] Dragomir Z. Djokovic. Corrections for "The closure diagram for nilpotent orbits of the split real form 

of £:8" [Cent. Eur. J. Math. 1 (2003), no. 4, 573-643 (electronic).]. Gent. Eur. J. Math., 3(3):578-579 

(electronic), 2005. 



18 W.A. DE GRAAF, E.B. VINBERG, AND O.S. YAKIMOVA 

[16] L.Yu. Galitski and D.A. Timashev. On classification of metabelian Lie algebras. J. of Lie Theory^ 

9:125-156, 1999. 
[17] The GAP Group. GAP - Groups, Algorithms, and Programming, Version 4-4: 2004. 

(http : //www . gap-system . org) . 
[18] Willcm A. de Graaf. Constructing algebraic groups from their Lie algebras, J. of Symb. Gomp., 44:1223- 

1233, 2009. 
[19] Willcm A. dc Graaf. SLA - computing with Simple Lie Algebras, a GAP package, 2009. 

(http: //science .unitn. it/~degraaf /sla.html) . 
[20] Willem A. de Graaf. Computing representatives of nilpotent orbits of 0-groups. J. of Symb. Gomp., 

46:438-458, 2011. 
[21] Fritz Grunewald and Joyce O'Halloran. Varieties of nilpotent Lie algebras of dimension less than six. J. 

Algebra, 112(2):315-325, 1988. 
[22] Sigurdur Helgason. Differential geometry, Lie groups, and symmetric spaces, volume 80 of Pure and 

Applied Mathematics. Academic Press Inc. [Harcourt Brace Jovanovich Publishers], New York, 1978. 
[23] V.G. Kac. Some remarks on nilpotent orbits. J. Algebra, 64(1): 190-213, 1980. 
[24] V.G. Kac. Automorphisms of finite order of semisimple Lie algebras. Funkcional. Anal, i Prilozen., 

3(3):94-96, 1969. 
[25] G. Kempken. Eine Darstellung des Kochers Ak. Dissertation, Rheinischen Friedrich-Wilhelms- 

Universitat Bonn, 1981; Bonner Mathematische Schriften 137(1982). 
[26] P. Littelmann. An effective method to classify nilpotent orbits. In Algorithms in algebraic geometry and 

applications (Santander, 1994), Progr. Math., 143, 255-269, Birkhauser, Basel, 1996. 
[27] Kenzo Mizuno. The conjugate classes of unipotcnt elements of the Chevalley groups Ej and £'§. Tokyo 

J. Math., 3(2):391-461, 1980. 
[28] D.D. Pervushin, Invariants and orbits of the standard (SL4(C)xSL4(C)xSL2(C))-module. Izv. Ross. 

Akad. Nauk Ser. Mat., 64(5):133-146, 2000; English translation in: Izv. Math., 64(5):1003-1015, 2000. 
[29] D.D. Pervushin, On the closures of the nilpotent orbits of fourth-order matrix pencils. Izv. Ross. Akad. 

Nauk Ser. Mat., 66(5):183-192, 2002; English translation in: Izv. Math., 66(5):1047-1055, 2002. 
[30] V.L. Popov. Two Orbits: When is one in the closure of the other? Tr. Mat. Inst. Steklova 264(2009), 

Mnogomernaya Algebraicheskaya Geometriya, 152-164; English translation: Proc. Steklov Inst. Math. 

264(2009), no. 1, 146-158. 
[31] D. M. Snow. Weyl group orbits. ACM Trans. Math. Software, 16(1):94-108, 1990. 
[32] Nicolas Spaltenstein. Classes unipotentes et sous-groupes de Borel, volume 946 of Lecture Notes in 

Mathematics. Springer- Verlag, Berlin, 1982. 
[33] M. Vergne, Cohomologie des algebres de Lie nilpotentes. Application a I'ctude de la variete des algebres 

de Lie nilpotentes. Bull. Soc. Math. France, (98):81-116, 1970. 
[34] E.B. Vinberg. The Weyl group of a graded Lie algebra. Izv. Akad. Nauk SSSR Ser. Mat, 40(3):488-526, 

1976. English translation: Math. USSR-Izv. 10, 463-495 (1976). 
[35] E.B. Vinberg. Classification of homogeneous nilpotent elements of a semisimple graded Lie algebra. 

Trudy Sem. Vektor. Tenzor. Anal., (19):155-177, 1979. English translation: Selecta Math. Sov. 6, 15-35 

(1987). 
[36] E.B. Vinberg, V.V. Gorbatsevich, and A.L. Onishchik. Lie Groups and Lie Algebras, III. Sovrem. probl. 

matematiki. Fundam. napravL, t.41. Moskva: VINITI, 1990. English translation: Encyclopaedia Math. 

Sci.; vol. 41, Springer, Berlin, 1994. 
[37] E.B. Vinberg and V.L. Popov. Invariant theory. In Algebraic geometry, IV. Itogi Nauki i Tekhniki, 

pages 137-314. Akad. Nauk SSSR Vsesoyuz. Inst. Nauchn. i Tekhn. Inform., Moscow, 1989. English 

translation: Encyclopaedia Math. Sci., vol. 55, Springer, Berlin, 1994. 
[38] E.B. Vinberg and A.G. Elasvili. A classification of the three-vectors of nine-dimensional space. Trudy 

Sem. Vektor. Tenzor. Anal, 18:197-233, 1978. English translation: Selecta Math. Sov., 7, 63-98, (1988). 
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Dimension 
80 

79 

78 

77 

76 

75 

74 

73 

72 

71 

70 

69 

68 

67 

66 

65 

64 

63 

62 

61 



5 6 




65 67 



66 68 



60 69 70 

Figure 1 . Hasse diagram of nilpotent orbits in the case of trivectors; top half 
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Dimension 
62 

61 

60 

59 



64 63 



65 67 




58 



57 



56 



55 



54 



53 



52 



51 



49 



48 

45 
42 
38 
37 
36 

35 

30 100 

19 101 

Figure 2. Hasse diagram of nilpotent orbits in the case of trivectors; bottom half 
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Table 4: Characteristics of the nilpotent orbits in the 
case of 3-vectors. 



no. 


characteristic 


no. 


characterisitc 


1 


66 66 6 66 12 


2 


66606666 


3 


66606066 


4 


60606606 


5 


06060606 


6 


6 15 15615 


7 


60606006 


8 


06006060 


9 


00600606 


10 


6 15 15 15 


11 


15 114 15 1 


12 


22222422 


13 


15 15 15 


14 


22222222 


15 


114 115 14 


16 


60006006 


17 


30330603 


18 


06000600 


19 


22022222 


20 


6 105015 


21 


03030330 


22 


14011411 


23 


14 1113 11 


24 


20420604 


25 


00600006 


26 


6 10 14 105 


27 


30303030 


28 


040 20420 


29 


112 11114 


30 


15 00105 


31 


22220202 


32 


10 5 10 14 


33 


20220204 


34 


00006000 


35 


114 10 10 4 


36 


20202022 


37 


00 10500 1 


38 


2 1111112 


39 


10 14010 


40 


10 104011 


41 


30330003 


42 


30003030 


43 


110 13 10 1 


44 


03000303 


45 


10 112 111 


46 


3 102120 


47 


20420 004 


48 


11111111 


49 


20004020 


50 


02002202 


51 


110 112 11 


52 


2 10 3 110 


53 


11110 112 


54 


03003000 


55 


00030300 


56 


20200202 


57 


00000006 


58 


04002000 


59 


00300030 


60 


111110 11 


61 


10 20301 


62 


0001005 


63 


1110 1021 


64 


12 10 20 


65 


12 110 11 


66 


02002002 


67 


10 110 3 10 


68 


100014 


69 


110 110 11 


70 


100 1004 


71 


02200020 


72 


2 10 110 1 


73 


10 110 110 


74 


100 100 13 


75 


03000003 


76 


120 00 102 


77 


00200400 


78 


10 10 12 


79 


02000004 


80 


110 10 10 1 


81 


00200200 


82 


110 00103 


83 


00030000 


84 


10 10 10 2 
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Characteristics. 


85 


102000 1 


86 


2 10 10 02 


87 


0002000 2 


88 


10 10 11 


89 


00000030 


90 


200 10 100 


91 


200 10003 


92 


10 100101 


93 


000 10020 


94 


1000 110 


95 


100 100 10 


96 


00000200 


97 


00100100 


98 


30000000 


99 


20000010 


100 


10001000 


101 


00 100 00 







Table 5: Characteristics of the nilpotent orbits. 



no. 


characteristic 


no. 


characterisitc 


1 


10 10 10 10 10 10 10 20 


2 


10 10 10 10 10 10 10 


3 


10 10 10 10 10 10 


4 


10 10 10 10 10 


5 


10 10 10 10 


6 


1 9 1 9 10 10 1 9 


7 


10 0377373 


8 


3 7373433 


9 


10 10 10 


10 


00 100 00 10 


11 


1 9 1 10 1 9 


12 


20822262 


13 


03703073 


14 


0000000 10 


15 


1 8 1 1 10 1 1 8 


16 


3 343 33 13 


17 


50500550 


18 


12 713 16 3 


19 


05050505 


20 


10 0019 


21 


22242422 


22 


60400640 


23 


33133124 


24 


03074303 


25 


10 10 118 


26 


10 10 10 9 


27 


0820 10 208 


28 


42224024 


29 


08200020 


30 


22222222 


31 


12 164312 


32 


110 110 18 


33 


10 111117 


34 


32053502 


35 


32052030 


36 


11111116 


37 


20000208 


38 


02002026 


39 


00500500 


40 


3 12 13 12 1 


41 


22064402 


42 


12 13 2 112 


43 


2 10 110 8 


44 


3304 1030 


45 


05000005 


46 


0003300 7 


47 


1111112 4 


48 


110 12 116 


49 


02222202 


50 


04000006 


51 


20202024 


52 


00300304 


53 


1112 1114 


54 


13 110 113 


55 


1023016 


56 


12 112 04 


57 


12 110 114 


58 


03200203 


59 


111112 13 


60 


10 113 106 


61 


02022022 


62 


00000050 
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63 


2 112 1111 


64 


02200204 


65 


11112 112 


66 


2 12 1112 


67 


00 100 140 


68 


0300003 1 


69 


50005000 


70 


10 110 4 


71 


12 10 12 1 


72 


00050000 


73 


2 10 13 12 


74 


10 10 1130 


75 


30201020 


76 


11111111 


77 


02000040 


78 


10400 1 


79 


00204006 


80 


2 110 10 2 1 


81 


110 10 13 


82 


00040002 


83 


10 11112 


84 


1202011 


85 


1030011 


86 


20002030 


87 


20111101 


88 


00202020 


89 


02 00 02 10 


90 


2 10 10 3 


91 


10 12 10 1 


92 


00030020 


93 


110 11110 


94 


10 2 110 


95 


1000300 


96 


10 1110 10 


97 


10 112 


98 


3000000 1 


99 


00020200 


100 


200 100 10 


101 


002000 10 


102 


10 1110 


103 


10 10 10 


104 


000 12000 


105 


100 1000 







Table 6: Characteristics of the nilpotent orbits. 



no. 


characteristic 


no. 


characterisitc 


1 


6 6 6 6 6 6 12 


2 


6660666 


3 


060 6 66 6 


4 


06 06 6 12 


5 


242 646 6 


6 


6060606 


7 


6 15 15 15 


8 


0606006 


9 


060 66 


10 


40 242 6 6 


11 


006 00 6 


12 


15 15 6 1 


13 


15 11451 


14 


22 2 2442 


15 


2040 20 6 


16 


2222222 


17 


150115 


18 


242 046 


19 


0600066 


20 


20 2 22 24 


21 


114 1114 


22 


2 13 1115 


23 


0006060 


24 


0303006 


25 


220 2 22 2 


26 


10 23133 


27 


0000600 


28 


3033003 


29 


0600000 


30 


10506 1 


31 


12 10 233 


32 


12 3 115 


33 


2042 004 


34 


1401111 


35 


0010501 


36 


0303030 
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37 


202 22 2 


38 


040 20 2 


39 


10 14 10 


40 


22 0242 


41 


0000006 


42 


0200400 


43 


2 111112 


44 


02 040 24 


45 


10 10 4 11 


46 


0202022 


47 


110 13 1 


48 


000 10 15 


49 


10 112 11 


50 


3000330 


51 


10 10 14 


52 


20 042 


53 


1020130 


54 


0300003 


55 


3010220 


56 


0200202 


57 


2010310 


58 


02 00 04 


59 


002 22 


60 


10 0203 


61 


10 024 


62 


110 1111 


63 


10 1022 


64 


0003000 


65 


1001112 


66 


10 2001 


67 


10 110 10 


68 


0002002 


69 


10 10 11 


70 


200 1000 


71 


0020000 


72 


10 10001 


73 


0000030 


74 


000 10 20 


75 


100010 
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Dimension 

48 

47 
46 
45 
44 
43 
42 

41 



40 38 



39 47 



6 7^ 10 /9 



11' 13. 12. 14 



15 17. 18._.19.^20 :i6 




38 

37 
36 
35 

34 

33 / 76' 7^ n77 73\ ^75 ^78(4,5) 

32 79 5i 80 82' 

Figure 3. Hasse diagram of nilpotent orbits; top 
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33 

32 79 



31 
30 
29 

28 
27 
26 
25 
24 
22 
21 



14 



76^ 74 77 73 .75 
80' 82' 




78 (4,5) 



95 (5,2) 



97 (5,2) 



104 (4,1) 



11 105(2,1) 

Figure 4. Hasse diagram of nilpotent orbits; bottom 
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Dimension 

42 

41 
40 
39 




38 



37 



36 



35 



34 



33 



32 41 



31 



30 



29 



28 



27 61 



26 



25 

23 

22 

21 69 

Figure 5. Hasse diagram of nilpotent orbits, E-j, top 
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Dimension 
29 



55^53 52 




11 75 (2,1) 

Figure 6. Hasse diagram of nilpotent orbits, Ej, bottom 



